Growth dynamics of a Bose-Einstein condensate in a dimple trap without cooling 
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We study the formation of a Bose-Einstein condensate in a cigar-shaped three-dimensional harmonic trap, 
induced by the controlled addition of an attractive "dimple" potential along the weak axis. In this manner we 
are able to induce condensation without cooling due to a localized increase in the phase space density. We 
perform a quantitative analysis of the thermodynamic transformation in both the sudden and adiabatic regimes 
for a range of dimple widths and depths. We find good agreement with equilibrium calculations based on self- 
consistent semiclassical Hartree-Fock theory describing the condensate and thermal cloud. We observe there is 
an optimal dimple depth that results in a maximum in the condensate fraction. We also study the non-equilibrium 
dynamics of condensate formation in the sudden turn-on regime, finding good agreement for the observed time 
dependence of the condensate fraction with calculations based on quantum kinetic theory. 
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I. INTRODUCTION 

The formation of a Bose-Einstein condensate (BEC) and 
the growth of long-range coherence from a gas of ther- 
mal atoms is a problem of interest in the field of ultracold 
atoms 0]]. Before the first observations of Bose-Einstein con- 
densation in a dilute gas there was some disagreement about 
the expected time scale for condensate formation J2l- Th e 
first quantitative predictions were made by Gardiner et al. JH, 
who derived a rate equation for the growth of a single con- 
densate mode from a super-critical thermal vapor. This was 
soon followed by the first experimental measurements of con- 
densate formation by Miesner et al. 0]. Starting from just 
above the critical temperature for a BEC, this experiment im- 
plemented a sudden evaporative cooling ramp to remove the 
high-energy tail of a near-degenerate Bose gas. The ensuing 
rethermalization led to the formation of a Bose-Einstein con- 
densate. These experiments were subsequently analyzed us- 
ing improved formalisms by Gardiner and co-workers JH-Q] 
and Bijlsma et al. QSQ . They found that, while their numerical 
calculations were qualitatively in agreement with experimen- 
tal observations, quantitatively no agreement could be found, 
and this has remained unresolved. 

In 2002 Kohl et al. [9] performed an experiment similar to 
that of Miesner et al. |4|] but implemented continuous rather 
than sudden evaporation from near quantum degeneracy. For 
this experiment the data were generally in good agreement 
with the results of quantum kinetic calculations incorporat- 
ing the details of the evaporation and the effects of three-body 
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loss lllOll . The same formalism was applied to later experi- 
ments in a quasi-condensate geometry by Hugbart et al. lUlll . 
where the calculated shape of the condensate growth curves 
agreed well with experiment apart from an unexplained time 
delay. Other evaporative cooling experiments leading to BECs 
worth noting are those of Schvarchuck et al. 11211 . who per- 
formed shock cooling in an elongated geometry and observed 
nonequilibrium dynamics in the resulting quasicondensate, 
and Ritter et al. IU3I1 . who measured the dynamics of the onset 
of long-range coherence in a three-dimensional condensate. 

Bose-Einstein condensation without evaporative cooling 
was first induced by Stamper- Kurn et al. l\M . motivated by 
the earlier work of Pinkse et al. Il5ll . Stamper- Kurn et al. 
began with a near-degenerate Bose gas in a cigar-shaped har- 
monic trap and slowly applied an additional attractive "dim- 
ple" trap formed by a red-detuned optical dipole potential to 
adiabatically increase the phase-space density by a factor of 
50. It was shown that this was reversible within the limits of 
heating caused by their dipole trap. Condensation was also in- 
duced by distillation without cooling in a double-well poten- 
tial, demonstrated in an experiment by Shin et al. Ilcjll . The 
lowering of a second well in this system caused the conden- 
sate in the first well to evaporate and re-form in the second at 
a higher temperature. Erhard et al. observed the formation of 
an rriF = BEC through spin collisions in an F = 1 spinor 
condensate from initial partically condensed components in 
the rriF = ±1 states, and they modeled their experiment us- 
ing rate equations fl7ll . Recently an alternate approach to re- 
versible BEC formation was demonstrated in an experiment 
by Catani et al. flUl . where entropy was exchanged between 
two atomic species, instead of between atoms inside and out- 
side a dimple potential. 

In this paper we revisit the method of Stamper-Kurn et 
al. ITill to quantitatively study the dynamics and the thermo- 
dynamics of Bose-Einstein condensation, and we make com- 
parisons of our experimental results with theoretical calcula- 
tions. We induce condensate formation by the controlled ap- 
plication of a tightly focused laser sheet to a near-degenerate 
Bose gas in a cigar-shaped magnetic trap (illustrated in Fig.Q]). 
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The addition of the resulting one-dimensional dimple poten- 
tial to the weakly confined dimension of the harmonic trap 
induces condensation by locally increasing density while the 
temperature remains almost constant, hence increasing the lo- 
cal phase-space density. 

We divide our results into two sections. First, we have mea- 
sured the final equilibrium state of the Bose gas following both 
the quasistatic (i.e., slow) and sudden turn-on of the dimple 
potential for a range of laser intensities beginning from a well- 
controlled initial nondegenerate state. The thermodynamics 
for an ideal gas with a delta function dimple have previously 
been studied in Ref. Il9ll . Using a self-consistent mean-field 
model involving semiclassical Hartree-Fock theory for the 
thermal cloud and the Thomas-Fermi approximation for the 
condensate, we can predict the final condensate fraction for a 
given dimple depth for both the quasistatic and sudden turn- 
on. We perform a quantitative comparison of experiment and 
theory for a thermodynamic transformation through the BEC 
phase transition in an interacting Bose gas. Second, we ob- 
serve the dynamics of condensate formation following sud- 
den turn-on of the dimple potential and compare with a quan- 
tum kinetic model of condensate growth. This configuration 
allows a quantitative comparison with theory for condensate 
formation without evaporative cooling. We note that this sce- 
nario has been studied previously using stochastic classical 
fields in one dimension Il20ll2lll . 

This paper is organized as follows: In Sec. [II] we summa- 
rize our experimental setup and procedure. In Sec. [Ill] we 
present our study of equilibrium thermodynamics by compar- 
ing the theoretical predictions with our experimental data and 
discussing the results. In Sec. IIV1 we present our study of 
condensate formation dynamics, first providing details of our 
quantum kinetic theory, and then comparing the theoretical 
predictions with our experimental data and discussing results. 
We finish with conclusions in Sec. [V] and we provide addi- 
tional theoretical details in the appendices. 



II. EXPERIMENTAL PROCEDURE 

Our experiments are performed on a gas of ultracold 87 Rb 
atoms in the h 2 Si/2\F = l,mj? = —1} state, confined in a 
magnetic harmonic potential provided by a Z-wire configura- 
tion on an atom chip, as previously described in Ref. 12211 . 
Briefly, we collect ~ 10 7 atoms in a reflection MOT be- 
fore transferring them to a magnetic trap at a distance of 
200 /im below the chip surface with axial and transverse 
frequencies of our cylindrically symmetric trapping potential 
(uj z ,u>±) « 2ir x (20,1600) Hz. Using rf evaporation we 
cool to near quantum degeneracy before decreasing the bias 
magnetic field to move the trap to 430 /im below the chip 
surface, with a resulting decrease in trapping frequencies to 
(uj z ,oj±) fa 2tt x (6.8, 160) Hz such that the resulting sys- 
tem is more three dimensional and less susceptible to phase 
fluctuations that exist in elongated condensates in the tighter 
trap 1I23I1 . Further evaporative cooling results in a cloud of 
(~ 2 - 6) x 10 5 atoms. We go as close as we can to the 
BEC transition while remaining above the transition temper- 




FIG. 1. (Color online) Schematic of the experiment. We create a 
dimple potential by tightly focusing a laser sheet at the center of the 
cigar-shaped harmonic magnetic potential. The total external poten- 
tial V C xt(z) along the axial dimension then consists of a wide har- 
monic potential plus a narrow Gaussian dimple potential, while the 
potential in the other two spatial dimensions is unaffected. 



ature. This ensures the largest possible condensate fraction 
when the dimple potential is subsequently turned on. 

Starting from these initial conditions, we next apply a red- 
detuned A = 840 nm optical dipole potential to the system, 
intersecting the weak (z) direction of the magnetic trap at the 
center as illustrated in Fig.[T] This dipole beam is known by 
direct charge-coupled device (CCD) imaging to be Gaussian 
and diffraction limited. By observing the effect of shifting 
the focus back and forth along the y dimension, we ensured 
the magnetic trap intersects the beam waist. The beam can 
be focused to two different 1/e 2 half-widths in the z dimen- 
sion, thereby creating either a wide (32 /j,m) or a narrow (11 
/im) dimple potential. This should be compared to the typical 
thermal cloud extent of ~ 400 /im near the critical tempera- 
ture. A cylindrical lens is used to expand the beam to widths 
of 350 /im (wide) or 220 /im (narrow) in the perpendicular 
dimension compared to the cloud width of 6 /im, such that the 
intensity is approximately constant in these dimensions. Thus 
the cylindrically symmetric trapping potential can be approx- 
imated by 



V cxt (r,z,t) = ^m(uj 2 j_r 2 + lu 2 z z 2 ) 



A(t)e 



-2{z/ w y 



(l) 



where w is the 1/e 2 half-width and the optical potential depth 



A( t) = i Io(t ) y ; klk - 



(2) 



is proportional to the peak laser intensity Io(t), and can reach 
a maximum depth of A{t)/k-Q = 210 nK (wide) or 1610 
nK (narrow). The other relevant parameters are the scatter- 
ing cross sections 07., linewidths 7^, the resonant frequencies 
u>k of the D\ and D2 lines in 87 Rb, and the laser frequency 
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Our measurements are performed using absorption imaging 
after a time of flight of 20.3 ms, after turning off all trapping 
potentials. The temperature is determined by fitting a thermal 
cloud distribution to the wings of the resulting image, and the 
condensate fraction is determined from a two-component fit 
to the density. The condensate is sufficiently dense that con- 
densate fractions of less that 1 % can be distinguished from the 
thermal background. 



III. THERMODYNAMIC TRANSFORMATIONS ACROSS 
THE CRITICAL POINT 

In our first set of experiments, we apply dimples of various 
depths and measure the resulting final equilibrium tempera- 
tures and condensate fractions. We do this for both wide and 
narrow dimples, and for both quasistatic and sudden turn-on. 
We then compare our measured values with the predictions of 
semiclassical Hartree-Fock theory, incorporating full mean- 
field interactions of both the condensate and thermal cloud, as 
well as accounting for the effects of three-body loss. While 
it is widely assumed that this is appropriate for the three- 
dimensional Bose gas, there have been relatively few com- 
prehensive comparisons with experimental data aside from 
Ref. & 



A. Theoretical procedure 

We determine the initial entropy Si and total energy E\ prior 
to dimple turn-on [A(t = 0) = 0], given the experimen- 
tally measured initial temperature and atom number, using the 
semiclassical theory as outlined in Appendix [A] To predict 
the final state, in the case of quasistatic turn-on we assume 
the system evolves isentropically and solve for the final tem- 
perature at which 5 = 5;. In the case of sudden turn-on, we 
use the initial density to calculate the sudden change in energy 
imparted by the dimple, 
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AE 



drAe- 2 ^ 2 [n (r) + n th (r)}, 



(3) 



where no(r) and nth(r) are the densities of the condensate 
and thermal cloud, respectively. We then solve for the final 
temperature at which E = E± + AE, which determines the fi- 
nal condensate fraction. We also estimate the effects of three- 
body loss in our calculations, as detailed in Appendix IB1 



B. Comparison with experiment — wide dimple 

Our experimental procedure is as follows. For the wide 
dimple measurements, we begin with an atomic cloud of N = 
6.25(25) x 10 5 atoms at a temperature of T x = 215(2) nK. 
This corresponds to a phase-space density at the center of the 
trap of r « 2.6, indicating that the cloud is very close to the 
BEC transition (r = C(3/2) ~ 2.612). To turn on the dim- 
ple quasistatically, we ramp up the dimple potential linearly at 
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FIG. 2. Comparison of theory and experiment for (a) the final tem- 
perature and (b) the condensate fraction at final equilibrium after the 
turn-on of the wide dimple potential. Quasistatic turn-on is indi- 
cated by solid lines for theory and by diamonds for experimental 
data. Sudden turn-on is indicated by dashed lines for theory and 
by squares for experimental data. All experimental data points are 
four-measurement averages with the error bars indicating the stan- 
dard deviation of the mean. 



a rate of 70/cb nK/s, and then hold it constant for a 300-ms 
equilibration time before turning off all potentials and imag- 
ing. To ensure that we are in the quasistatic regime, we have 
repeated this process for various ramp rates: For faster rates 
we observe a decrease in condensate fraction at large dimple 
depths due to nonadiabatic heating. To turn on the dimple 
suddenly, we ramp up the potential in less than 0. 1 ms, and 
then hold it constant for a 1000-ms equilibration time before 
imaging. 

The resulting temperatures and condensate fractions are 
plotted versus final dimple depth and compared with the pre- 
dictions of semiclassical theory in Fig. 12 We find good agree- 
ment between theory and experiment for both the quasistatic 
and sudden turn-on, with the best fits obtained using initial 
conditions N = 6.50 x 10 5 and 1] = 216 nK. We also find the 
results of the quasistatic and sudden turn-on to be nearly in- 
distinguishable from each other. This is because only a small 
fraction of the atoms are drawn from the harmonic trap into 
the wide dimple, which has a maximum attainable depth of 
order ~ fee J 1 . 

We compare the predictions of the semiclassical theory 
with our experimental data for both wide (Fig. 12) and narrow 
(Fig.0) dimples. We plot temperature and condensate fraction 
versus dimple depth for both quasistatic and sudden turn-on. 
In all cases, the only fitting parameters used to generate the 
theoretical curves are the initial temperature and atom number, 
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FIG. 3. Comparison of theory and experiment for (a) the final tem- 
perature and (b) the condensate fraction at final equilibrium after the 
turn-on of the narrow dimple potential. Quasistatic turn-on is in- 
dicated by solid lines for theory and by diamonds for experimental 
data. Sudden turn-on is indicated by dashed lines for theory and by 
squares for experimental data. The dotted lines indicate quasistatic 
turn-on from the same initial conditions used for sudden turn-on. All 
experimental data points are single measurements, with the exception 
of a four-point measurement at depth A/ks = 1610 nK to measure 
the standard deviation of the mean. 



sistatic turn-on, with the best fits obtained using initial condi- 
tions N = 2.45 x 10 5 and 7] = 162 nK. 

However, we do not find good agreement for the sudden 
turn-on, with theory predicting a much smaller condensate 
fraction than was experimentally observed. The best fit shown 
is for N — 2.25 x 10 5 and T\ = 166 nK, though better fits can 
be obtained by using values of N and 2] that lie outside their 
respective measurement uncertainties. 

A potential explanation for this discrepancy is that the turn- 
on is not sufficiently quick, and the density of the gas is not 
frozen during the turn-on. This would result in a smaller in- 
crease in total energy of the gas compared to the prediction 
of our model, leading to a smaller increase in temperature 
and hence a larger condensate fraction — as is observed ex- 
perimentally. Indeed, in the limit of quasistatic turn-on from 
the same initial conditions (dotted lines in Fig. [3j, the pre- 
dicted condensate fraction is considerably larger than for sud- 
den turn-on. Since the measured values lie between the lim- 
its of quasistatic and sudden turn-on, an intermediate turn-on 
time would likely provide a good fit to the data. 

To model this would require a fully dynamical treatment 
that is beyond the limitations of our semiclassical theory and 
is numerically impractical within the quantum kinetic model 
we introduce later in Sec. llV Al Furthermore, departures from 
the sudden case would only be expected at turn-on times sim- 
ilar to the time scale for rethermalization. Since the mean free 
time between collisions at initial conditions is approximately 
6 ms, whereas the time for sudden turn-on is less than 0.1 
ms, this explanation seems unlikely. We are therefore unfor- 
tunately forced to leave this discrepancy unresolved. 



constrained to lie within their respective measurement uncer- 
tainties. 



D. Discussion 



C. Comparison with experiment — narrow dimple 

Using the more tightly focused narrow dimple, we are 
able to attain a maximum depth much greater than k^Ti and 
thereby observe differences between quasistatic and sudden 
turn-on. For the quasistatic turn-on we begin with an atomic 
cloud of N = 2.60(15) x 10 5 atoms at a temperature of 
Ti = 160(2) nK (r w 2.6). We ramp up the dimple potential 
linearly over a time of 1500 ms, and then hold the potential 
constant for a 300-ms equilibration time before imaging. For 
the sudden turn-on we begin with an atomic cloud of N = 
2.10(15) x 10 5 atoms at a temperature of 2] = 168(2) nK 
(T « 1.1) and follow the same procedure as for the wide dim- 
ple. 

The resulting temperatures and condensate fractions are 
plotted versus final dimple depth and compared with theory 
in Fig. [3] In this case we observe a significant difference be- 
tween the results of quasistatic and sudden turn-on, though 
this is partly due to the difference in initial conditions. We find 
good agreement between theory and experiment for the qua- 



An interesting feature of the narrow dimple data is that there 
exists an optimal dimple depth (~ 750/cb nK) at which a max- 
imum condensate fraction is obtained. As the dimple depth is 
increased beyond this value, the condensate fraction gradually 
decreases back toward zero. We can understand this feature in 
the context of the condensate formation process as follows. At 
shallow depths, the dimple potential acts merely as a perturba- 
tion to the broader harmonic trap. To a first approximation the 
chemical potential and temperature are unchanged, whereas 
the energy of the translational ground state is decreased rela- 
tive to the bottom of the harmonic trap (see Fig.|4|. When the 
ground-state energy approaches the chemical potential, a con- 
densate forms, as observed for the wide dimple and at shallow 
depths of the narrow dimple. However, at larger depths a sig- 
nificant fraction of the thermal cloud is drawn into the narrow 
dimple, causing non-negligible changes in /x and T. At suf- 
ficiently large depths, the entire thermal cloud falls into the 
dimple potential (which is approximately harmonic near its 
center) and the condensate evaporates, in agreement with the 
well-known result that compression of a gas cannot alter the 
phase-space density ll26ll . 
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FIG. 4. Schematic of the effect of dimple turn-on for small depths A. 
Before the dimple is turned on (solid lines), as the system tempera- 
ture is above T c the chemical potential n lies below the translational 
ground state energy eo of the harmonic trap. When the dimple is 
turned on (dashed lines), the new ground state energy e' lies below 
the initial chemical potential, inducing condensation. Due to interac- 
tions, the energy e' increases as the condensate grows until equilib- 
rium is reached (e' « fi). 



IV. CONDENSATE FORMATION DYNAMICS 

In our second set of experiments, we suddenly turn on the 
dimple and allow the system to evolve for various times before 
measuring the condensate fraction. We do this for both wide 
and narrow dimples, in each case for two different fixed dim- 
ple depths: one shallow and the other deep. We then compare 
our measured values with the predictions of quantum kinetic 
theory. In particular, we use the ergodic quantum Boltzmann 
equation 12711 with the additional inclusion of the effects of the 
condensate mean-field J3, [H an d three-body loss. While sim- 
ilar comparisons have been made previously 10-Qj]] the con- 
densation transitions in these experiments were induced by 
evaporative cooling. Here condensation is introduced without 
loss by modifying the density of states of the trap. 



A. Quantum kinetic theory 

The initial equilibrium state of the Bose gas above T c be- 
comes nonequilibrium following the sudden turn-on of a dim- 
ple potential. We begin by calculating this initial nonequilib- 
rium state, and then evolve it in time to final equilibrium using 
an ergodic quantum Boltzmann equation (EQBE), as previ- 
ously described in 

The ergodic assumption is that all semiclassical states for a 
given single-particle energy e have the same mean occupation 
and thus the nonequilibrium phase-space distribution function 
depends only on time and energy: /(r, p, t) = /(e(r, p, t),t). 
Hence the position and momentum dependence in the full 
quantum Boltzmann equation is projected out, yielding a par- 
tial differential equation in energy and time only. The ergodic 
approximation is necessary in order to reduce the dimension- 
ality of the phase space to render the quantum Boltzmann 
equation computationally tractable. The EQBE then governs 



the evolution of the energy density 

n(£ ' t)= I (2^^ £ ~ £th(r ' P ' <} ) /(r ' P ' t] (4) 
= g(e,t)f(e,t), 

where f(e, t) is the energy distribution function and 

5(M)= / (2^ <5(£_eth(r ' P '* )) (5) 

is the density of states. The semiclassical excitation en- 
ergy eth(r, p, t) is as defined in Eq. dAU but with the mean 
field of the thermal cloud neglected [i.e., setting n t h(r) to 
zero]. This eliminates the need to determine the densities self- 
consistently, and our results suggest that this approximation is 
reasonable. The condensate mean field is calculated via the 
Thomas-Fermi approximation, as per Eq. ( |A3t , but again ne- 
glecting the mean field of the thermal cloud. We provide a 
more detailed description of the EQBE in Appendix [C] and 
describe how to incorporate the effects of three-body loss in 
Appendix IB1 

A typical atom will travel a distance of less than 1 /im dur- 
ing the sudden ramp-up of the dimple potential. As this is 
more than an order of magnitude smaller than the width of 
the narrow dimple, we approximate the sudden turn-on as in- 
stantaneous. Under this assumption, the phase-space distribu- 
tion immediately after the sudden turn-on must be the same 
as before the turn-on: /'(r,p) = /(r,p) = / B E(e t h(r, p)), 
where we henceforth use primed (unprimed) variables to de- 
note quantities immediately after (before) dimple turn-on at 
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FIG. 5. Energy density n(e, t) = g(e, t)f(e, t) before dimple turn- 
on at t < (solid line), immediately after dimple turn-on at t = 
(dashed line), and at final equilibrium after t = 1000 ms (dotted 
line). Strictly speaking, the value plotted at e = — A is the conden- 
sate occupation, not the energy density. Immediately after dimple 
turn-on to depth A, only a small fraction of atoms have energies be- 
low the bottom of the harmonic potential (e < 0). Over the 1000-ms 
equilibration time, this fraction steadily increases, eventually result- 
ing in the formation of a condensate. 
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Time since dimple turn-on (ms) 



FIG. 6. Condensate growth curves for theory and experiment follow- 
ing the sudden turn-on of the wide dimple potential. The condensate 
fraction is plotted vs time for dimple depths of A/ks — 70 nK (solid 
line: theory; diamonds: experiment) and A/ks = 140 nK (dashed 
line: theory; squares: experiment). All experimental data points are 
four-measurement averages with the error bars indicating the stan- 
dard deviation of the mean. 

t = 0. Although the phase-space distribution function is un- 
changed, the post-dimple semiclassical excitation energy dif- 
fers from the pre-dimple expression via the inclusion of the 
dimple potential: 

e' th (r,p) = e th (r,p)-Ae- 2 ^ 2 . (6) 

Clearly, the change in e[ h (r, p) alters the density of states in 
Eq. ©, and therefore it alters the energy density, which we 
calculate immediately after the dimple turn-on as 

n'{e)= J (2^3 5 ( e ~ e tfr( r > P)) /be ( e th(r, P)) ■ (7) 

In doing this we are ergodically projecting the pre-dimple 
phase-space distribution using the post-dimple density of 
states. It should be noted that the physical initial phase space 
will actually be non-ergodic. However, previous Monte Carlo 
calculations of the Boltzmann equation have shown that er- 
godicity is restored relatively quickly J28ll . and hence we ex- 
pect this should be a reasonable approximation. 

The pre- and post-dimple energy densities are shown in 
Fig-El for parameters typical of the narrow dimple: A/k-g, = 
805 nK, N = 1.30 x 10 5 , and T; = 133 nK. The post- 
dimple energy density closely resembles the pre-dimple en- 
ergy density, except that the peak in post-dimple energy den- 
sity is slightly smaller because a small fraction of atoms lying 
within the dimple now have energies below the bottom of the 
harmonic potential (e < 0). Also plotted in Fig.|5]is the energy 
density at final equilibrium, calculated via the EQBE. As the 
system evolves toward equilibrium, successively more atoms 
are drawn into the dimple, occupying the energies in the range 
— A < e < 0, and ultimately resulting in the formation of a 
condensate. 
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FIG. 7. Condensate growth curves for theory and experiment follow- 
ing the sudden turn-on of the narrow dimple potential. The conden- 
sate fraction is plotted vs time for dimple depths of A/ks = 660 nK 
(solid line: theory; diamonds: experiment). All experimental data 
points are four-measurement averages with the error bars indicating 
the standard deviation of the mean. 



B. Comparison with experiment 

We compare the predictions of the quantum kinetic cal- 
culations just described with our experimental data for both 
wide (Pig. [6]) and narrow (Fig. [7]) dimples. In both cases, we 
plot condensate fraction versus time for two different dim- 
ple depths. In all cases, the only fitting parameters used to 
generate the theoretical curves are the initial temperature and 
atom number, constrained to lie within their respective mea- 
surement uncertainties. 

For the wide dimple measurements, we begin with an 
atomic cloud of N = 6.25(25) x 10 5 atoms at a temperature 
of T ; = 215(2) nK (r « 2.6). We ramp up the potential sud- 
denly, and then hold it constant for various equilibration times 
up to 1000 ms before imaging. We use dimples of two differ- 
ent depths: A/kg = 70 and 140 nK. The resulting condensate 
fractions are plotted versus time and compared with the pre- 
dictions of quantum kinetic theory in Fig. [6] We find good 
agreement between theory and experiment, with the best fits 
obtained using initial conditions (N = 6.30 x 10 5 , 2] = 217 
nK) and (N = 6.10 x 10 5 , T- = 216 nK) for dimple depths 
A/k-Q = 70 and 140 nK, respectively. 

For the narrow dimple measurements, we begin with an 
atomic cloud of N — 1.87(7) x 10 5 atoms at a temperature 
of T; = 154(6) nK (r « 1.4), and we follow the same exper- 
imental procedure described earlier for a final dimple depth 
of A/k-Q = 660 nK. The results are plotted in Fig. |7l with 
the best fit obtained using initial conditions N — 1.94 x 10 5 
and Tj = 148 nK. We find good agreement between theory 
and experiment, particularly at final equilibrium. In contrast 
to the wide dimple case, we observe a delay of about 50 ms 
in the onset of condensate formation, which agrees with the 
EQBE calculation. 
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V. CONCLUSIONS 

We have studied the thermodynamics of the transition 
across the BEC critical point and the formation dynamics of 
Bose-Einstein condensation resulting from the application of 
attractive Gaussian dimple potentials of various widths and 
depths to a Bose gas in a cigar-shaped harmonic trap. We 
have measured the equilibrium temperature and condensate 
fraction over a range of dimple depths, for both quasistatic 
and sudden turn-on of the dimple, and compared our results 
with the predictions of semiclassical Hartree-Fock theory. For 
narrow dimples, we found that there exists an optimal dimple 
depth at which a maximum condensate fraction is attained. 
Beyond this depth, the dimple acts merely as a tighter har- 
monic trap, and it therefore does not increase the phase-space 
density. We found good agreement between theory and exper- 
iment, except in the case of sudden turn-on of deep, narrow 
dimples. We also measured the (nonequilibrium) condensate 
fraction over a range of times after sudden turn-on of the dim- 
ple and compared our results with the predictions of quantum 
kinetic theory. We observed a short delay in the onset of con- 
densate formation in the case of narrow dimples, but not in 
the case of wide dimples. In both cases we found good agree- 
ment between theory and experiment. This provides further 
validation to the quantum kinetic model of condensate forma- 
tion, for which previous comparisons were based on sudden 
evaporative cooling. 
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a. The mean-field density combined with the external poten- 
tial Vgxt (r) constitute the effective potential experienced by 
atoms in the thermal cloud. The density of the thermal cloud 
is calculated as 



nth(r) 



dp 



(2nh) 



3/BE(eth(r,p)), (A2) 



while the density of the condensate is calculated using the 
Thomas-Fermi approximation, 



(A3) 



ra (r) = maxjo, (/x - [V cx t(r) + 2U n t h(r)] 

Thus, the densities of both the thermal cloud and conden- 
sate must be solved self-consistently for a given T and fx to 
give the experimentally measured total atom number N = 
J dr [n (r) + n t h(r)]. 

From here other thermodynamics quantities such as the 
condensate fraction Nq /N, total entropy, and total energy can 
be determined. The last two are relevant to quasistatic and 
sudden turn-on of the dimple potential, respectively. By writ- 
ing^ e = e t h(r, p) the total entropy is given by (cf. p. 15 of 



drdp 



/BE(e) + llln[/BE(e) + l 



-/BE(e)ln /be(c) 



while the total energy of the system is 

E = K + V + 1, 
drdp p 2 



K = 



-/be(c), 



(27rft) 3 2m' 
V= I drF ext (r)[no(r)+nth(r)], 



(A4) 

(A5) 
(A6) 

(A7) 



I=^-Jdr [n 2 (r) + 4n (r)n th (r) + 2n? h (r)] .(A8) 



Appendix A: semiclassical theory 

We make use of the Thomas-Fermi approximation for the 
condensate density and the semiclassical Hartree-Fock ap- 
proximation for the thermal cloud, as outlined in I29T43511 and 
compared with experiment by Gerbier et al. II25I1 . 

At equilibrium the thermal cloud is well described 
by the Bose-Einstein distribution, /BEUthC 1 ', p)) = 
(cxp{[e t h(r, p)— n]/k-B,T} — 1) _1 , where T is the temperature 
and p, is the chemical potential. In the semiclassical Hartree- 
Fock approximation, the excitation energy e t h(r, p) (the en- 
ergy required to add an atom with momentum p at position r) 
is given by the expression 

P 2 

e th (r, p = i- + Fext(r) + U [2n (r) + 2n th (r)] , (Al) 
2m 

where m is the atomic mass and Uq — Anh 2 a/m is the inter- 
action strength, proportional to the s-wave scattering length 



Appendix B: Three-body loss 

Three-body loss occurring between condensed and non- 
condensed atoms can be calculated from the three-body cor- 
relation function as 

n(r) 3B = -/v 3 (^ t (r) 3 *(r) 3 ), (Bl) 

where K 3 = 5.8(1.9) x 1CT 30 cm 6 /s J37t]. Using a bro- 
ken symmetry approach, we write the Bose field operator as a 
mean field plus fluctuations, 

£(r) = V(r) + <5(r), (B2) 

and substitute this into the previous expression. Identifying 
the condensate density as no(r) = |f/>(r)| 2 yields 

n(r) 3B = -K 3 {[n (r)} 3 + 9[n (r)} 2 {S\r)S(r)) 

+9n (v)(S\r) 2 6(v) 2 ) + (5t(r) 3 5(r) 3 )}.(B3) 
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The noncondensate density is given by nth(r) = (^( r )^( r ))- 
Using Wick's theorem on the higher order operator moments 
of the fluctuations gives 



(*t(r) 2 i(r) 2 )=2[ )lth (r)] 2 , (B4) 

(Jt( r ) 3 5(r) 3 ) =6[n th (r)] 3 , (B5) 

-A' 3 {Mr)] 3 +9Mr)] 2 n th (r) 

48n (r)[n th (r)] 2 + 6[n th (r)] 3 }. (B6) 



and so 



h{r) 



.'SB 



The loss rates for the condensate and thermal cloud atoms can 
then be written separately by noting that, for example, the sec- 
ond term in this above expression represents a three-body pro- 
cess in which two condensate atoms and one thermal cloud 
atom are lost. The coefficients can be divided accordingly to 
giveJH 

™o(r) 3B = -A 3 {[7i (r)] 3 +6[n (r)] 2 n th (r) 

+6n (r)[n th (r)] 2 }, (B7) 
nth(r) 3B = -if 3 {3[n (r)] 2 n th (r) + 12n (r)[n th (r)] 2 

+6[n th (r)] 3 }, (B8) 

which ensures no(r) 3 B + n-th(r)3B = «(r) 3 B- To incorporate 
these loss rates in our semiclassical theory (Appendix lAl. we 
treat these expressions for rto( r )3B and rith(r) 3 B as coupled 
ordinary differential equations, and we solve them in time at 
each spatial gridpoint via Euler's method. This is not entirely 
straightforward, because the semiclassical theory is strictly 
static: The final equilibrium is determined from the initial 
equilibrium in a single leap, without any stepwise time evolu- 
tion. In the case of sudden dimple turn-on, we begin with the 
final equilibrium densities, and we calculate three-body loss 
over the 1-s equilibration time. In essence, we are assuming 
that the three-body loss rates on the final equilibrium densi- 
ties are approximately equal to the average loss rates during 
the 1-s evolution from initial to final equilibrium. In the case 
of the quasistatic dimple turn-on, we can be more precise: Be- 
cause the system is never out of equilibrium as it evolves from 
initial to final equilibrium, we can break up the process into 
arbitrarily many steps. For each step, we evolve the system 
isentropically as the dimple depth is incrementally increased, 
and we calculate three-body loss over the corresponding time 
interval. 

The inclusion of three-body loss has a significant effect on 
the predictions of both semiclassical Hartree-Fock and quan- 
tum kinetic theory, as we show in Fig. [8] This is particularly 
evident in the case of adiabatic turn-on, where the inclusion 
of three-body loss not only drastically reduces the maximum 
condensate fraction but also reduces the optimal dimple depth 
at which the maximum fraction is attained. All of the curves 
shown correspond to the narrow dimple, with initial condi- 
tions N = 2 x 10 5 and T ; = 150 nK. 
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FIG. 8. The effects of three-body loss on the predicted condensate 
fraction, (a) Semiclassical Hartree-Fock predictions of the equilib- 
rium condensate fraction vs dimple depth in a narrow dimple for adi- 
abatic turn-on, including three-body loss (solid line), adiabatic turn- 
on, with no three-body loss (dotted line), sudden turn-on, including 
three-body loss (dashed line), and sudden turn-on, with no three- 
body loss (dot-dash line), (b) Quantum kinetic theory predictions of 
condensate growth following the sudden turn-on of the narrow dim- 
ple of depth A/ks = 660 nK, including three-body loss (solid line) 
and with no there-body loss (dashed line). 



Appendix C: Ergodic quantum Boltzmann equation 

In our implementation of the EQBE, we use an energy grid 
with uniform spacing 5e, with the lowest energy bin corre- 
sponding to the condensate. Because the condensate energy 
level (equal to the bottom of the effective potential) changes 
with time, we redefine our energy grid to be e' = e — (i[No(t)], 
where n[No(t)] is the condensate chemical potential in the 
Thomas-Fermi approximation. With this transformation the 
EQBE is d] 



h(e,t) 



d_ 
"de 



ffw(e,*) 



n(e,t) 



+ I co ii(e,t), 



where <7 w (e, t) is the weighted density of states, defined as 



(CI) 



drdP ^( £ -e th (r,p,i))^M, (C2) 



Of 



which contains the time derivative of the effective potential 
V(r, t) = K, x t(r) + 2[/rj7io(r, t) and hence depends on the 
rate of condensate growth N (t) = h(0,t) = l co \\(0,t). The 
contribution of binary elastic collisions is given by the term 



m 3 g 2 



27r 3^7 / de2de 3 de4,g(€mm)S(£ + £2 - £3 - £4) 

X [(1 + /)(1 + /,:/:,/; - // 2 (1 + / 3 )(1 + / 4 )] , 

(C3) 



where / = /(e,t), fi = f(ei,t), and e min = 
min{e, 62, e 3 , £4}. The first (second) term within the square 
brackets represents the forward (backward) collisions e 3 + 
£4 <-> £ + £2- Factors of the form (1 + fi) are due to Bose 
enhancement and vanish in the classical limit /; <C 1. For 
the case of collisions involving a condensate atom, / C oii(0, t), 
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we make the replacement g(0)[l + f(0,t)} w g(0)f(0,t) = 
N (t). Detailed derivations of the EQBE can be found 
in SHH. 

We numerically evolve the EQBE in time using an explicit 
fourth-order Runge-Kutta method. At each time step, we cal- 
culate Icon by summing over all combinations of energy bins 
{e, €2, £3, £4} satisfying the delta function in Eq. (|C3t and up- 
dating I co ii(ei, t) for all four bins involved in each combina- 
tion. We then add the contributions from three-body loss [see 
Eq. ( |C6t and dC7H , and use the resulting value of -ZV (£) to 
calculate g w (e, t) in the first term on the right-hand side of the 
EQBE. 

To incorporate three-body loss in our quantum kinetic cal- 
culations, the loss rate for the thermal cloud must be modified 
to include energy dependence: 

«th(e,r) 3B = -i^ 3 n th (e,r)|3[n, (r)] 2 

+ 12no(r)nth(r) + 6[n th (r)] 2 }, (C4) 

where the energy-position density of the thermal cloud is cal- 
culated from the energy density as 



n th (e,r) = g(e,r)f(e) = ^-n(e), 



and g(e, r) is the position-dependent density of states, ob- 
tained by omitting the spatial integral from Eq. Q. At each 
time step in the numerical evolution of the EQBE, we deter- 
mine the condensate density no(r) from the Thomas-Fermi 
approximation [Eq. iA3i . with the thermal cloud mean-field 
neglected] and determine the thermal cloud density n t h(r) by 
integrating this expression for n t h(e, r) over energy. We then 
calculate the loss rate for the condensate, Eq. ( lB7t . and the 
energy-dependent loss rate for the thermal cloud, Eq. iC4\ , 
Lastly, we integrate out the spatial dependence to get 



(C5) 



N 03 b = J drh (r) 3B , (C6) 
"th(e)3B = J drn t h(e,r)3B, (C7) 



which are simply added to the EQBE calculation of n(e) in 
Eq. CD- 
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